Generic Gravitational Wave Signals from the Collapse of Rotating Stellar Cores 



O 

o 

(N 



(N 
> 

in 
o 

m 

(N ■ 

o ; 

i> 

o 

Oh 

I ' 

o ■ 



H. Dimmelmeier/'^ C. D. Ott,^ H.-T. Janka,^ A. Marek,^ and E. Miiller^ 

^Department of Physics, Aristotle University of Thessaloniki, GR-54124 Thessaloniki, Greece 
^ Max-Planck-Institut fur Astrophysik, Karl-Schwarzschild-Str. 1, D-85741 Garching, Germany 
'^Department of Astronomy and Steward Observatory, University of Arizona, Tucson, AZ 85721, USA 

We perform general relativistic (GR) simulations of stellar core collapse to a proto-neutron star, 
using a microphysical equation of state (EoS) and an approximation of deleptonization. We show 
that for a wide range of rotation rates and profiles the gravitational wave (GW) burst signals from the 
core bounce are generic, known as Type I. In our systematic study, using both GR and Newtonian 
gravity, we identify and quantify the influence of rotation, the EoS, and deleptonization on this 
result. Such a generic type of signal templates will facilitate a more efficient search in current and 
future GW detectors of both interferometric and resonant type. 

PACS numbers: 04.25.Dm, 04.30.Db, 95.30.Sf, 95.30.Tg, 95.55.Ym, 97.60.Bw 



Theoretical predictions of the gravitational wave (GW) 
signal produced by the collapse of a rotating stellar iron 
core to a proto-neutron star (PNS) in a core collapse 
supernova are complicated, as the emission mechanisms 
are very diverse. While the prospective GW burst signal 
from the collapse, bounce, and early postbounce phase 
is present only when the core rotates 0, 0, S, 0, 01' 
one expects GW signals with sizeable amplitudes also 
from convective motions at later post-bounce phases, 
anisotropic neutrino emission, excitation of various os- 
cillations in the PNS, or nonaxisymmetric rotational in- 
stabilities BBBEIill- 

In the observational search for GWs from merging bi- 
nary black holes, powerful data analysis algorithms like 
matched filtering are applied, as the waveform from the 
inspiral phase can be modeled very accurately. In stark 
contrast, the GW burst signal from stellar core collapse 
and bounce cannot yet be predicted with the desired 
accuracy and robustness. First, a general relativistic 
(GR) description of gravity and hydrodynamics including 
the important microphysics is necessary. Only very few 
multi-dimensional codes have recently begun to approach 
these requirements. Second, the rotation rate and profile 
of the progenitor core are not very strongly constrained. 
Therefore, the influence of rotation on the collapse dy- 
namics and thus the GW burst signal must be investi- 
gated by computationally expensive parameter studies. 

Previous simulations, considering a large variety of ro- 
tation rates and profiles in the progenitor core but ignor- 
ing complex (though essential) microphysics and/or the 
infiuence of GR, found qualitatively and quantitatively 
diff'erent types of GW burst signals (see, e.g., 0, 0, S])- 
These can be classified depending on the collapse dynam- 
ics: Type I signals are emitted when the collapse of the 
homologously contracting inner core is not strongly influ- 
enced by rotation, but stopped by a pressure-dominated 
bounce due to EoS stiffening at nuclear density pnuc 
where the adiabatic index 700s rises above 4/3. This leads 
to instantaneous PNS formation with a maximum core 
density Pmax > Pnuc- Type II signals occur when centrifu- 



gal forces grow sufficient to halt the collapse, resulting in 
consecutive (typically multiple) centrifugal bounces with 
intermediate coherent re-expansion of the inner core, seen 
as density drops by often more than an order of mag- 
nitude; thus here Pmax < Pnuc after bounce. Type III 
signals appear in a pressure-dominated bounce when the 
inner core has a very small mass at bounce due to a soft 
subnuclear EoS or very efficient electron capture. 

In contrast, new GR simulations of rotational core col- 
lapse employing a microphysical EoS and an approxima- 
tion for deleptonization during collapse suggest that the 
GW burst signature is exclusively of Type I f/J . In this 
work we considerably extend the number of models and 
comprehensively explore a wide parameter space of ini- 
tial rotation states. Also for this more general setup we 
find GW signals solely of Type I form. We identify the 
physical conditions that lead to the emergence of this 
generic GW signal type and quantify their relative influ- 
ence. Our results strongly suggest that the waveform of 
the GW burst signal from the collapse of rotating super- 
nova cores is much more generic than previously thought. 

We perform all simulations in 2 -|- 1 GR using the 
CoCoNuT code approximating GR by the con- 

formal flatness condition (CFG), whose excellent quality 
in the context of rotational stellar core collapse has been 
demonstrated extensively (see, e.g., 0, EH)- CoCoNuT 
utilizes spherical coordinates with the grid setup speci- 
fied in [7] and assumes axisymmetry. GR hydrodynam- 
ics is implemented via finite- volume methods, piecewise 
parabolic reconstruction, and an approximate Riemann 
solver. We extract GWs using a variant of the Newtonian 
quadrupole formula (see, e.g., [l3|). 



We employ the microphysical EoS of Shen et al. [14 



15| . Deleptonization by electron capture onto nuclei and 
free protons is realized as in [l^: During collapse the 
electron fraction Yg is parameterized as a function of den- 
sity based on data from spherically symmetric ncutrino- 
hydrodynamic simulations [l5| using the latest available 
electron capture rates 17 1 (updating recent results Q 
where standard capture rates were used). After core 
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FIG. 1: Time evolution of the GW amplitude h and maximum 
density pmax for three representative models with different 
rotation profiles and initial rotation rates /3i. 



bounce, Y^. is only passively advected and further lepton 
loss is neglected. Contributions due to neutrino pressure 
Pi, are taken into account above trapping density [3| ■ 

As initial data we take the non-rotating 20 Mq progen- 
itor s20 from imposing the rotation law discussed 
in 0, @] ■ In order to determine the influence of different 
angular momentum on the collapse dynamics, we param- 
eterize the initial rotation of our models in terms of the 
differential rotation parameter A (Al: A ~ 50,000 km, 
almost uniform; A2: A = 1,000 km, moderately differ- 
ential; A3: A = 500 km, strongly differential) and the 
initial rotation rate /3i — T/\W\, which is the ratio of 
rotational energy to gravitational energy (approximately 
logarithmically spaced in 18 steps from 0.05% to 4%). 

As in Q , in the entire investigated parameter space our 
models yield GW burst signals of Type I, i.e. the wave- 
form exhibits a positive pre-bounce rise and then a large 
negative peak, followed by a ring-down (upper panel of 
Fig. [T]) . However, with respect to collapse dynamics and 
the relevant forces halting the collapse, the models fall 
into two classes. While for instance all models with the 
almost uniform rotation profile Al experience a pressure- 
dominated bounce for which a Type I waveform is ex- 
pected, models with profiles A2 or A3 and sufficiently 
strong initial rotation rate (3i (> 4% for A2; > 1.8% 
for A3) show a single centrifugal bounce at subnuclear 
density (lower panel of Fig. [T]). Nevertheless, they also 
produce a Type I waveform, as their core does not re- 
expand after bounce to densities much less than bounce 
density but immediately settles to a PNS after a short 
ring-down phase. What obviously distinguishes models 
with pressure-dominated bounce from those with cen- 
trifugal bounce is that the latter have GW signals with 
significantly lower average frequencies. Note also that 



models with very little rotation develop convective over- 
turn of the shock-heated layer immediately after shock 
stagnation (not to be confused with the late-time con- 
vection discussed in Q), resulting in a lower-frequency 
contribution to the post-bounce GW signal (see Fig[T|). 

In order to analyze the absence of Type II signals, 
in particular for cases with centrifugal bounce, we now 
separately investigate the influence of GR, a microphys- 
ical EoS, and deleptonization on the dynamics of rota- 
tional core collapse with different amounts and distribu- 
tions of angular momentum. When the pre-coUapse iron 
core starts to contract, its effective adiabatic index 7ofr 
is lower than the critical value ~ 4/3 needed for stabil- 
ity against gravitational collapse. Here jeS is the sum of 
the adiabatic index 7cos = dlnP/d\np\Y^^s of the EoS 
(where P is the pressure, p the density, and s the spe- 
cific entropy of the fluid) and a possible correction due to 
deleptonization (which can be significant until neutrino 
trapping; see 0). At this stage, both GR and rotational 
effects (in our range of /3i) are negligible in discussing 
stability. If the build-up of centrifugal forces in the in- 
creasingly faster spinning core during collapse is strong 
enough, contraction is halted and the core undergoes a 
centrifugal bounce rather than reaching nuclear density 
(where EoS stiffening would also stop the collapse). 

A necessary condition for a centrifugal bounce at sub- 
nuclear densities is that "f^s exceeds a critical rotation 
index 7rot- There exists a simple Newtonian analytic re- 
lation, 7rot = (4 - 10/3ic,b)/(3 - 6/?ic,b) [H (where /3ic,b 
is the inner core's rotation rate at bounce), which works 
well in equilibrium, but is rather imprecise as a crite- 
rion for centrifugal bounce in a dynamical situation. For 
instance, for rotating core collapse models in Newtonian 
gravity with a simple hybrid EoS [So^ and no deleptoniza- 
tion (where 7cff = 7cos), we find that for our range of ini- 
tial rotation rates and 1.24 < 7cff < 1.332, the analytic 
relation strongly underestimates the actual 7rot by up to 
~ 0.2 at high /3ic,b- Furthermore, /3ic,b is a result of the 
evolution, depending on the initial parameters A and /3i 
of the prc-coUapse core in an a-priori unknown way. 

For this reason, we determine the boundary between 
pressure-dominated and centrifugal bounce in terms of 
the parameters (3i and 7cff by systematic numerical sim- 
ulations. For models with a simple hybrid EoS ^20] and 
no deleptonization, using the same initial density profile 
as in the microphysical models, the results are shown in 
Fig. [21 both in the Newtonian case (dashed lines) and 
in GR (solid lines). Apparently, for our choice of ini- 
tial rotation the influence of GR can be approximated 
by adding an offset of — A7gi. ~ 0.015 to the Newtonian 
results (dotted lines). This gives a quantitative measure 
of the GR effects on rotational core collapse, which is in 
agreement with [4]. Note that A7gi- is negative because 
GR effectively acts like a softening of the EoS. 

Fig. [2] also shows for each rotation profile the loca- 
tions where the transition between pressure-supported 
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FIG. 2: Boundary between pressure-dominated and centrifu- 
gal bounce in the 7eff-/3i plane for models using the hy- 
brid EoS in Newtonian gravity (dashed lines) and GR (solid 
lines). Dotted lines show the Newtonian results shifted by 
— A7gr — 0.015. The transition points for models using the 
microphysical EoS without and with deleptonization, again 
for Newtonian gravity (circles) and GR (bullets), lie in the 
shaded areas around 7off — 1.32 and 1.29, respectively. 



and centrifugal bounce occurs when the microphysical 
EoS is used and Pi is gradually increased. These tran- 
sitions are marked on the different boundary lines and 
allow the identification of the 7eff value where models 
with the hybrid EoS make this transition. For simula- 
tions with microphysical EoS but no deleptonization we 
find that in all cases the transition occurs near 7eff ~ 1.32 
(highlighted by the upper grey band in Fig.[2]) . This value 
agrees with the average of 7eos for the microphysical EoS 
at densities between 10^^ and 10^^ g cm~^, which is the 
most relevant range for the collapse dynamics. Thus the 
type of bounce obtained with the microphysical EoS is 
well reproduced by a hybrid EoS with 7cos — 1.32. 

Deleptonization before neutrino trapping reduces 
7off compared to 7eos locally according to Aj^ — 
^S\nYe/S\np < [21[. Above trapping an additional 
positive correction Aj^ « 5{P^ / P) / 6h\ p due to neutrino 
pressure effects must be considered. From the Yf.-p tra- 
jectories used to describe the deleptonization during core 
collapse, we estimate values between —0.03 and —0.02 
for A7e -l- again in the density regime relevant for 
the bounce dynamics. Therefore 7cff « 1.29 is expected 
for models with microphysical EoS and deleptonization. 
Again this agrees with the simulations. Fig. [2] shows that 
the circles and bullets marking those models on the differ- 
ent boundary lines (for the investigated initial rotation 
profiles with Newtonian gravity and GR) all lie in the 
range of values indicated by the lower grey band. 

The finding that deleptonization decreases 7cff to 
about 1.29 explains the absence of Type II GW signals 
for our most sophisticated models. When a hybrid EoS 
is used, the subgroup of such models showing multiple 
centrifugal bounces and subsequent strong re-expansion 
phases of the inner core occupies only a small area in the 



7cff-/3i plane. It is located at 7eff > 1-31 for all of our 
initial rotation states both in the Newtonian case and in 
GR, i.e. significantly above the value of 7cff ~ 1.29 that 
characterizes the microphysical models if deleptonization 
is included. 

The value 7off = 1.29 captures well the deleptonization 
effects in the density regime above neutrino trapping, and 
therefore models with hybrid EoS and 700s — 1-29 exhibit 
the same collapse and bounce behavior in this dynamical 
phase. For 7cos = 1.29 the hybrid EoS corresponds to 
an effective average Yg of « 0.237, and thus for models 
with this EoS we find a small mass Afjc ^ 0.1- 0.3 Mq 
of the inner core at bounce (higher for more rapid ro- 
tation), consistent with the theory of self-similar col- 
lapse |22j . For models with microphysics, however, 7off is 
much closer to 4/3 in the early collapse phase where the 
initial value of Afjc is determined. At intermediate den- 
sities below and around neutrino trapping, deleptoniza- 
tion indeed reduces Mic, but only to about 0.5 - 0.9 M©, 
which is in agreement with recent spherically symmet- 
ric GR results using Boltzmann neutrino transport 
These large values of Mic explain the complete absence 
of Type III GW burst signals in our models (see also [^). 

The generic nature of the GW burst signal has several 
important implications for prospective detectability. The 
limitation to a unique Type I waveform for a very broad 
range of rotation states of the progenitor core will very 
likely facilitate the use of more powerful and finetuned 
data analysis methods in GW detectors. To this end, 
we offer our results in a waveform catalog 2j|. Note 



that almost all investigated models result in a pressure- 
dominated bounce and instantaneous formation of a PNS 
with similar average density and compactness. For these 
models, whose rotation rates at bounce span two orders 
of magnitude (0.2% < /3ic,b ^ 20%), the maxima /max of 
their waveform's frequency spectrum lie in a very narrow 
range with an average of /max — 718 Hz. 

This clustering in frequency is reflected in Fig. [31 where 
we plot the (detector-dependent) frequency-integrated 
characteristic waveform amplitude hr. against the char- 
acteristic frequency /c (Eq. (31) in [25|) for all 54 GR 
models with microphysical EoS and deleptonization. We 
assume optimal orientation of source and detector, and in 
cases with pressure-dominated bounce remove the lower- 
frequency contribution from the post-bounce convective 
overturn by cutting the spectrum below 250 Hz, as we 
are only interested in the GW signal from the bounce 
and ring-down. Fig. [3] shows that while current LIGO 
class interferometric detectors are only sensitive to sig- 
nals coming from an event in the Milky Way, advanced 
LIGO could marginally detect some signals from other 
galaxies in the Local Group like Andromeda. For the 
proposed EURO detector .^2Q\ in xylophone mode, we ex- 
pect a very high signal-to-noise ratio (which is he divided 
by the detector sensitivity at /c). This detector could also 
measure many of the computed signals at a distance of 
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FIG. 3: Location of the GW burst signals from the core 
bounce for all models in the /ic~/c plane relative to the sensi- 
tivity curves of various GW interferometer detectors (as color- 
coded). The sources are at a distance of 10 kpc for LIGO, 
0.8 Mpc for advanced LIGO, and 15 Mpc for EURO. Open 
circles denote models where the filtered-out early postbounce 
convection contributes significantly to the original signal, and 
crosses show models undergoing centrifugal bounce. 



15 Mpc, i.e. in the Virgo cluster. 

Note that detectabiUty could be enhanced by (i) a net- 
work of interferometers, (ii) the support by resonant de- 
tectors, which is facilitated by the narrow range of /max, 
and (iii) the use of more powerful data analysis meth- 
ods based on the waveforms' similarity and robustness. 
A serious obstacle for detection is the low event rate of 
~ 1 yr~^ within 15 Mpc, which is further reduced by as- 
suming that only a fraction of all progenitors rotate fast 
enough to have a strong GW bounce signal. Neverthe- 
less, our results can serve as a guideline for a possible fre- 
quency narrow-banding of future interferometers (which 
can significantly boost sensitivity), as well as for choos- 
ing the optimal configuration in the planning of resonant 
detectors. In addition, as a core collapse may be accom- 
panied by other GW emission mechanisms like late-time 
convection (also in a nonrotating core) or bar-mode insta- 
bilities, the total GW signal strength and duration could 
be significantly higher than predicted here. 

As a downside, the generic properties of the GW burst 
signal introduce a frequency degeneracy into the signal 
inversion problem. Consequently, in the case of a de- 
tection it is difficult to extract details about the rota- 
tion state of the pre-coUapse core, because a large part 
of the corresponding parameter space yields a pressure- 
dominated bounce and thus signals with very similar val- 
ues for /max and (depending on the detector) also for 
/c. On the other hand, as /max directly depends on the 



compressibility of the nuclear EoS at bounce, determin- 
ing this frequency from the GW burst signal can help to 
constrain the EoS properties around nuclear density. 
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